The Scalings of Scalar Structure Functions in a Velocity Field with Coherent Vortical 

Structures 
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In planar turbulence modelled as an isotropic and homogeneous collection of 2-D non-interacting 
compact vortices, the structure functions S p (r) of a statistically stationary passive scalar field have 
the following scaling behaviour in the limit where the Peclet number Pe — * oo 

for LPe~ 1/3 < r < L , 



<S p (r) ~ constant + In 



(LPe" 1 / 3 ) 



S P (r) 



LPe- 1 / 3 



6(1-D) 



for LPe~ 1/2 < r < LPe" 1/3 , 



where L is a large scale and D is the fractal co-dimension of the spiral scalar structures generated 
by the vortices (1/2 < D < 2/3). Note that LPe' 1 / 2 is the scalar Taylor microscale which stems 
naturally from our analytical treatment of the advection-diffusion equation. The essential ingredients 
of our theory are the locality of inter-scale transfer and Lundgren's time average assumption. A 
phenomenological theory explicitly based only on these two ingredients reproduces our results and a 
generalisation of this phenomenology to spatially smooth chaotic flows yields (fclnfc) -1 generalised 
power spectra for the advected scalar fields. 



I. INTRODUCTION 

The theory of turbulent passive scalars has received 
much attention recently I^Qjj- The mixing of a scalar 
held 9 in a velocity held v is governed by the advection- 
diffusion equation 

9 t 0(x, t) + v(x, t) ■ V0(x, t) = kV 2 6»(x, t) + /(x, t) (1.1) 

where k denotes the molecular diffusivity of the scalar 
9 and /(x, t) is an external source ( forcing) driving the 
scalar. The mixing of a scalar held 9 is characterised 
by its structure functions S p = ([9(x + r) — 9(x)] p ) = 
((59(r)) p ) for any number p. If we want to hnd S p then 
we need models of the velocity held v. 

The model which has attracted much attention re- 
cently is the Kraichnan model Q where the velocity 
held v is considered to be incompressible, statistically 
isotropic, white-noise in time(J-correlated) and Gaussian. 
Furthermore it has homogeneous increments with power 
law spatial correlations 

( [ Vi (r , t) - v t (0, 0)] [ Vj (r, t) - Vj (0, 0)]) = 



25{t)r h (h + d-l)Sij-h^ 



(1.2) 



where the scaling exponent h S]0, 2[ and d is the di- 
mension of space so that i,j = 1,2, ... ,d. The above 
tensorial structure of the velocity held is in conformity 
with incompressibility. The Kraichnan model also as- 



sumes a forcing /(x, t) in (1.1) that is an independent 
Gaussian random held with zero mean. The forcing is 
white in time and its covariance is assumed to be a real, 



smooth, positivc-dehnite function with rapid decay in 
space so that the forcing is homogeneous, isotropic and 
takes place on the integral scale L. The generic scaling 
behaviour of the structure function S p ~ , r <€. L, of 
passive scalars in the Kraichnan model was established in 
@'@0Jll • The scaling exponents of this formalism are of 
the form £ p = ( p (d, h,p) where h £ [0, 2] is the Holder ex- 
ponent in (1.2). In the context of this model Balkovsky 
and Lebedev || and Chertkov || used the instantonic 
formalism in a d- dimensional space to hnd the scaling 
exponents for large p. It was also shown in the instan- 
tonic formalism that linip^oo ( p ~ rf( - 2 gfe fe ' 1 [|| which is 
independent of p. The scaling exponents were also cal- 
culated using other techniques in the limits h — ► ||,^]], 
d — > oo [^]|7]] and p — > oo g] , and a 2 — h expansion of Cp 
was proposed in ||. It was found that £ p does depend 
on p for small values of p in the Kraichnan model. 

Our work lies in the opposite extreme of the Kraichnan 
model. We work in the regime where we have a persis- 
tent vortical velocity held frozen in time in two dimen- 
sions. The important differences between this model and 
the Kraichnan model are in the structure of the velocity 
held inhnitely correlated in time in this model but delta 
correlated in time in the Kraichnan model; and vortical 
in space in this model, but Gaussian in the Kraichnan 
model. In this model, the velocity held is that of spatially 
distributed nonintcracting two dimensional vortices with 
compact structure. We consider the spatial distribution 
of vortices to be dilute in that they are far from each other 
and therefore maintain their structure and spatial posi- 
tion for an indefinite period. We also consider this distri- 
bution to be homogeneous and isotropic and the velocity 
held to be incompressible, that is V • v = 0. The model 
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of the velocity field considered here is an artificial model 
of planar homogeneous turbulence where the emphasis is 
on the coherent vortex aspect of the flow. In order of pre- 
sentation, the first aim of this model is to demonstrate 



that in the case of the unforced scalar (/ = in 1.1) we 
can quantify the statistics of the turbulent scalar field 
in terms of the scalar's spiral geometry generated by the 
coherent vortical structures in the flow(sections |l| and 
IV). The second aim is to derive the Batchelor /c^^ower 
spectrum and all the corresponding structure functions 
for the scalar field in the case where the scalar is forced 
and in a statistically steady state (section |v|). Such a 
spectrum has been recently observed by Jullien et. al., 
|| in a 2-D turbulent flow with well defined, albeit short 
lived, coherent vortical structures. 

In the next section we discuss the scenario of a decaying 
scalar field in an isolated vortex. In section III , we define 



structure functions and calculate the spectrum of higher 
order correlation functions for a single spiral created by 
a single vortex and decaying by the action of molecu- 
lar diffusion. In section [[v] we generalise our analysis to 
many non-interacting vortices and calculate the scalings 
of the structure functions of the decaying scalar field. In 
section |v| we calculate the generalised power spectra and 
the corresponding structure functions of statistically sta- 
tionary scalar spirals by applying the time-average oper- 
ation approach of Lundgren . In section VI we discuss 
the phenomenology behind the fc _1 scalings of the gen- 



eralised power spectra. Section VII contains conclusions, 
discussion and the obtainment of the (fclnfc) -1 scalings 
of the generalised power spectra in smooth chaotic flows. 



II. PASSIVE SCALAR IN A PLANAR VORTEX 

The advection of a decaying passive scalar field by a 
single planar vortex has been studied by Flohr and Vas- 
silicos |ll[. We use the formulation used in namely: 



d t e + ft(0<V = kv 2 <? 



(2.1) 



d d - 

where dt = -rr and cL, = — ■ and Q(r) = fin (t) and 
dt d(j) VL/ 

L is the maximum distance of the scalar interface from 
the centre of the vortex. This equation describes the ad- 
vection and diffusion of a scalar field 9 in the azimuthal 
plane r = (r, <fi) by a steady vortex with azimuthal ve- 
locity component u<f,{r) = rd(r) — LUq (-j-) s and van- 
ishing radial and axial velocity components. Direct Nu- 
merical Simulations and experiments in the laboratory 
have demonstrated the existence and importance of co- 
herent vortices in two-dimensional turbulence and in two- 
component turbulence in stably-stratified flow with and 
without rotation of the reference frame ]l^-|l5|). Note 
that axial velocity fields of the form u^{r) — LSln (-£) 1 s 
have been used in Jl0|,O,|l6|-pTl and that their large 



wavenumber energy spectrum has the form E(k) ~ 
k~ 5+2s for 1/2 < s < 2 with the appropriate large scale 
bound. We choose s > 1 to ensure that u<p(r) does not 
increase with increasing r and s < 2 to ensure that the 
energy spectrum is steeper than k . The initial scalar 
field 6*o = #(x, t = 0) is characterised by a regular inter- 
face between 9q =/= and 9q = with minimal distance r$ 
and maximal L from the rotation axis. By regular struc- 
ture we mean that the interface has no irregularities on 
scales smaller than L. Nothing else needs to be specified 
about the initial scalar field 9q(x). Such a patchy initial 
condition where all the non-zero scalar is confined within 
a regular interface mimics well initial conditions in cer- 
tain laboratory experiments where scalar is injected in 
the flow in the form of blobs(e.g. ||). 

As time proceeds, the patch winds around the vortex 
and builds up a spiral structure and decays due to diffu- 
sion. The characteristic time Qq 1 is the inverse angular 
velocity of the vortex at L. This defines a Peclet nu mbe r 
Pe = VIqL 2 ^ 1 . We non-dimensionalise equation (2.1) 



by using the following transformations 



r, n a t -> t, x fi(r) fl(r), L 2 V 2 -> V 2 



and equation fl2.l|) takes the form 

1 



d t d + n{r)d^9 = j^-V 2 9 . 



(2.2) 



In the non dimensionalised variables we have Q(r) = r~ s , 
and ro represents rg/L since L = 1. Considering finite 
diffusivity k, the form of the general solution of the evo- 
lution equation ^. 2[ ) for any initial field 9q in the limit of 
large t, ie. t ^> 1, is the following [jll]] 

9(r 1 t) = J2fn(r,t)e m ^ n ^ 

n 

/„(r,*)=/ n (r,0)e[-i r 

(2.3) 

where r =\ r | and <fi is the azimuthal angle and n is an 
integer, il' is the derivative of fl with respect to r. The 
angular Fourier coefficients /„ (r, t) are time dependent 
and the initial condition is fully specified by / „(r, 0). We 
do not go in to the details of the solution of ( |2.l[) which 
can be found in |1]0|. 



III. STRUCTURE FUNCTIONS OF PASSIVE 
DECAYING SCALAR IN ONE VORTEX 

In this work we concentrate in finding the scaling prop- 
erties of the structure functions of the scalar field in a pla- 
nar turbulence consisting of coherent vortices. The two 
point equal time p-th order structure function is defined 
as follows 
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Sp(r,f) = ([0(x + r,t)-0(x,i)] p ) 

= PMF)- 



(3.1) 



Sp depends only on the magnitude of the distance be- 
tween two points, when the ensemble average is taken 
over an isotropic and homogeneous distribution of the 
scalar field. The overbar denotes ensemble averaging and 
the brackets imply space averaging ((...) oc f (fix). 

Let us first calculate ([S0(r, t)] p ) for one 2-D vortex. 
We use the binomial expansion as follows 

([0(x + r,t)-0(x,t)n 

= <^(x + r,t)) + (-l)^(x,t)) 

p-i 

+ £C' W (-l) 9 <0 P ~ , (* + r,t)0«(x,t)) 

q=l 

(3.2) 

where C qp is the binomial coefficient of the expansion. In 
order to calculate ([66(r, t)] p ) we first determine 

(0*(x, t)flP-«(x + r, *)) = B pq (r, t) (3.3) 

which is the g-th term in the binomial expansion of 



((S8(r)) p ) as shown in (3.2). Now we write the above 



as 



^ J d 2 x0«(x,t)0P-«(x + r,t)=5 M (r,*) 



(3.4) 



where La is a large scale over which the spatial average 
may be calculated. The Fourier transform of equation 
(3.4) is given by 



F pg Qt,t) 



1 

2^ 



-ik-r 



B pq (r,t)d 2 v. 



(3.5) 



Substituting (3.4) in (|3.5| ) and after some standard ma- 
nipulations we get 

4 9 (M) = ^r / 0«(x,t)e ik ' x d 2 x 



x / d p - q (x',t)e 



-ik-x' d 2 x , 



(3.6) 

Now if we integrate ( |3.6| ) over a circular shell in k space 
we get 



F pq {k,t) 



2tt 



dA(k)F pq (k, t) 



(3.7) 



where dA(k) = kd<pk, k — |k| and fa is the angle of k. 
F pq (k, t) could be called the generalised power spectrum 
of the scalar field in Fourier space. Substituting (2.3) and 
(3.6) in (3.7) we get the following 



F pq {k,t) 
1 



(2nL A ) 



dA(k) 



in(<l>-fl(x)t)] 



p-q 



,[im(0'-n(x')*)] 



(3.8) 



where x = |x| and x' = |x'|. After some standard ma- 
nipulations ( |3.8| ) leads to 

F Pg (k,t) 

/ / dxxJ n {kx)2Tr(i) n e m ' 1 "' 



(2nL\) 



^ ' fnifn-2 ■ ■ ■ fn—n-L. 



-in£l{x)t 



n,n 1 ,ri2...n q - 1 

x J dx'x'e^J^ikx'^ni-i)" 1 

x ^ ^ f mi f rri2 fm—mx...—m p -q-i 



x e 



m,mi,m2...m p _g_i 



(3.9) 



where we have changed summation variables in accor- 
dance with the following conditions 



n = Til + Tl2 + «3 • • • "■<; 

and to = toi + m 2 + to 3 + . . . TOp_ g . 



(3.10) 



All the /„'s and / m 's are functions of time i and of x and 
a/ respectively. J n (kx) is the Bessel function which has 
been substituted in place of the integral representation 



J,n4> -ikxcoa{4>-<t>k) J J, — Ott( i\ n e in( t>^ 



2Tr(i) n e m * k J n (kx). (3.11) 



After integrating ( |3.9| ) with respect to fa and summing 
over to we get the following relation 

F pq (k,t) 

= , n \ 9 s [ dxkxJ n (kx)2Tr(i) n 
{2nL A ) J 

s, f f f -inQ(x)t 

n,ni,n2...n ? _i 

X y dx'x',P n (kx')2TT(~i)- n 

x ^ ] /mi /m 2 • • • f-n-m 1 ...—m p _ q _ 1 



X e 



mi,m 2 ---m P -i3-l 



(3.12) 
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FIG. 1. The two cases (I & II) when the centre of the vortex 
is respectively outside a scalar patch and inside it. 

We now have two cases to study (see figure 1). In case 
I the vortex centre is outside the scalar patch and its 
nearest distance to the scalar interface is tq. In Case II 
the centre of the vortex is inside the scalar patch and its 
nearest distance to the scalar interface is tq. In view of 



the above we can write (3.12) as 



F pq (k,t) 



(2nL\) 



dxkxJ n [kx)2n(i) n 



^ ] /ni/«2 • • • fn— Tii...— n,_i ( 



-inQ,(x)t 



ro roo' 



J 7*0 



dx'x'J*_ n (kx')2ir{-iy 

^ ^ /mi/m.2 • ■ ■ f—n—mi...- 



x e 



mi,m 2 ...m p - g _i 



(3.13) 



In (3.13) we have four terms of the form as shown below 

dx'{...) 



dx 



ro rr 

dx / dx' (. . .) 
Jo 

/>oo pro />oo />oo 

+ dx dx'{...)+ / dx dx'{...) 

J ro Jo J ro J ro 

(3.14) 

with the integrands denoted by (...) being the same as 
( |3 .13 ) for all the terms of the above. For case I, terms 
containing Jq°(- ■ ■ ) are zero since /„ = in the regime 
< x < tq for all n, even for n — 0, since there is no 
scalar patch in the region < x < rQ. Therefore contri- 
butions only come from the term 



dxx / dx x (. . .). 

J To 



For case II we can legitimately replace 9 by 9 — fo and get 
the same result as in case I. Hence the only contributing 
term is the following 



F pq (k,t) 
1 



^L\> Jr a 



dxkxJ n (kx)2ir(i) n 



^ . fn 1 fn 2 ■ ■ ■ /n— Tii...— Tlg-i ( 



n,m,7i2...rig_i 



x / dx'x'Jl n {kx')2ir(-i)- n 

^ ^ frnifrri2 ■ ■ ' f—n—mi... — r 



X e 



m 1 ,m 2 ...m p _ g _i 
inQ(x')t 



(3.15) 



It is because / decays with a time scale which is much 
larger than the decay time of the non-zero harmonics |l^] 
that fo is considered to be a constant and therfore sub- 
tracted away from the 9 in case II. The same reasoning 
can be applied for case I. Hence all the n,'s and m,'s 
are non-zero in the above equation and in the rest of the 
paper . 

To take into account the fact that diffusion gradually 
wipes out the spiral structure of the scalar field near the 
vortex centre (a fact not taken in to account in |l6| ] where 
the spiral structure is assumed to exist wholly intact un- 
til finally destroyed by viscosity), we follow Flohr and 
Vassilicos [0 and define a critical radius p which gives 
a measure of this diffused region. This critical radius is 
defined in the limit Pe — ♦ oo for times t <C Pe 1 / 3 which 
are such that 



fn(r,t) = f„(r,0) forr<t:p 



but 



\fn{r,t)\ « |/„(r,0)| for r»p See (2.3) 



Hence, we set ^n 2 Vt l2 (p)Pe 1 t 3 = 1 which implies 



Pit) 



-n 2 s 2 Pe- 1 t 3 
. > 



This critical radius is time-dependent and grows with 
time. It can be thought of as a diffusive length scale over 
which the harmonics in n have diffused and the spiral 
structure has been smeared out. 



In view of the above, the integrals in ( 3.15 ) can be 
further divided as 



dx 



dx 



dx. 



The only significantly non-zero contribution comes from 
the range p < x < 1 in the integrals and similarly for x' . 
Hence we get the following 



F pq (k,t) 
1 



(?*L* A ) J p 

E 

n,ni,n 2 ...nq_i 



dxkxJ n (kx)27r(i) n 

fnifri2 • • • fn—n\...—n 



-in£l(x)t 
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x / dx'x'Jl n (kx')2n(~i)- n 
Jp 

^ /mi/m 2 • ■ ■ / — n— mi. 

mi,m2...m p -,-i 



(3.16) 



Now for large kx, ie. h> 1, we can use the asymptotic 
expansion for the Bessel function 



J n {kx) 



1 



2nkx 



(3.17) 



This is appropriate for our analysis if the Fourier modes 
are to resolve at the very least the distance ro from 
the centre of the vortex to the scalar patch interface 
i.e 1 < kr . After substituting ( 3 . 1 7| ) in ( 3.16 ) we use 
the method of stationary phase to evaluate the integrals 
where the phase is given by 



$ = kx — nfl(x)t. 



(3.18) 



The approximation for a general integral of this type is 
known to be 



I(x) = I f(t)exp[ix^(t)]dt 

-f{t*)exp[ix^{t*)±Tr/4] 



2x | |' 



(3.19) 



where t* is the t where the derivative of the phase is zero. 
The condition of stationary phase gives 



= = -k - nn'(x n )t 



(3.20) 



which picks out points x n where the contribution to the 
integral is maximum. Finally what we get is 



F pq {k,t) r 

E 

ro,m . ..n q -i 



27rfc 



n\W(x n )\t 



X V 2nk ) ^ ni ^ n2 ' ' ' ni ~ n2 ~ "' ~ n i- 1 



^ ^ fmi fm 2 



n—m\— rn 2 — ■ ■ ■ — Trip — q — i 



771 1 ,7712 ■ ■■Tn p _ q _i 



(3.21) 



where f n = f n (x n , 0) and simil arly for f m . Now from the 
condition of stationary phase ( 3.20 ) we can find 



' snt\ 3+1 



(3.22) 



where we have used f2(r) — r s . The stationary phase 
contributes only when p < x n < 1 because the spiral 



structure only exists in that range of distances from the 
centre of the vortex. The relation between the fractal co- 
dimension (Kolmogorov capacity) D of the scalar spiral 
and the power law of the decay of the azimuthal velocity 
of the vortex is El, p2| 



D = 



(3.23) 



This D is such that 1/2 < D < 2/3 because 1 < s < 2 
and gives a measure of the space-filling proper ty of the 
spiral. Hence after doing the summations in ( |3.21 ) we 
can show that the power spectrum F pq (k, t) scales like 

F pq (k,t) 

~ fc _ ( 3_2D ^ 2 ( 1_D ) [const + higher order terms] 

(3.24) 

in the limit Pe — * oo and in the range t < k < \J~^- 

for times 1 < t < Pe 1 / 3 which is the range of times for 
which the scalar patch has a well-defined spiral structure 

in the range of wavenumbers t < k < ^/^(obtained 
from p < x n < 1). The higher order terms are func- 
tions of k/t, and decay faster than (k/t) -1 in the range 

t < k < \J~^-, and can therefore be neglected. 

We notice that as time runs forward the spiral range 
(t < k < \p^p) shifts to higher values of k which is solely 
due to the vortex continuously wrapping the scalar field 
in to finner and finner spirals thus generating scales which 
have higher wavenumbers. This range also shrinks as it 
shifts to higher values of k because of the action of diffu- 
sion. In the next section we generalise our results to the 
case of multiple spirals generated by a dilute collection of 
nonintcracting vortices which may be representative of a 
turbulent velocity field with coherent vortical structure, 
perhaps obtained in the experiments of Jullien et al. [p|. 



IV. STRUCTURE FUNCTIONS OF PASSIVE 
DECAYING SCALAR IN A FLOW CONSISTING 
OF MANY IDENTICAL NON-INTERACTING 
VORTICES 

All the analysis in this section is carried out in dimen- 
sionalised variables so we invert the transformations of 
section 0. Let us consider many non interacting vortices 
randomly distributed over 2-D space and sufficiently far 
apart so that we can safely describe the velocity field in 
terms of compact vortical structures characterised by 



n(aO = n (f)" 



ft(x) 



if — > 1 

1j 
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where the x's are measured from the centre of each vortex 
at x m for all m and 



min\x m — x„| 3> L ; for all m and n. 



(4.1) 



For the calculation of the generalised power spectrum 
and structure functions we need only to consider the 
scalar patches within a distance L of each vortex be- 
cause these scalar patches acquire a spiral structure and 
thereby dominate the scaling of the structure functions. 
The scalar field at distances larger than L from all vor- 
tices contributes an 0(r/L) term to the structure func- 
tion for r -C L because the interfacial structure of the 
scalar field far from the vortices remain regular. As we 
show in this section, this term is negligible in the r/L<l 
limit compared to the r-dependence of the structure func- 
tions caused by the scalar spiral structures around the 
vortices. It is therefore sufficient to consider that 0(x, t) 
consists only of the local scalar fields # m (x — x m ,i) in 
the vicinity of vortices labelled m and write 



(x,t) = $^0 m (x- 



(4.2) 



Every scalar spiral in the right hand side of ( |4.2|) i s lo - 
calised within a distance L of x m and the condition (4.1) 
ensures that they do not ove rlap each other. 

Now we can generalise (3J3) to include the effect 
of many non-interacting vortices with non-overlapping 
scalar spirals as shown below 

4*(M) = ^E/ ^(x,t)e tkx d 2 x 

•A m 

x J ^-«(x',t)e-* x 'rf 2 x'. (4.3) 

Since the integrals are independent of the m's we can 
take them out of the sum. Hence we get the same result 
as in ( |3.24 ) multiplied by the number of vortices per unit 
area, that is 



F pq (k,t) ~ (kL) 



-(3-2D) 



(not) 



2(1-23) 



E 



i 



2nLY 



(4.4) 



This asymptotic relation is valid when Clot < kL < y 
which is found from the condition p < x n < I in di- 
mensionalised form and 1 <C flot <C Pe 1 / 3 in the limit 
Pe — > oo. 

Assuming the distribution of the scalar spirals over the 
two dimensional space to be homogeneous and isotropic, 
the power spectrum F pq (k,t) — 2?:kF pq (k,t) where the 
bar implies ensemble averaging over the distribution of 
many spirals (because vortices are non-interacting and 
spirals are therefore statistically independent from each 
other, it does make sense for the average over space al- 
ready included in the definition of F pq (k.,t) to be taken 
over an idealised space where there is only one spiral, and 



for the ensemble average to be taken over the distribu- 
tion of many spirals). From (3.2) we can show that for a 
homogenous distribution of scalar spirals the odd order 
structure functions vanish. Only the even order struc- 
ture functions do not vanish, that is for p — even. Hence 



from(3.1) 



S p (r,t) = ([0(x + r,t)-0(x,t)] 



= ({59(r,t)]P} 

-/(»-£<* 



■ik-r 



F pq (k, t)kdkd<j) 

(l-e lk - r )F pq (k,t)dkd(t) 
(1 - J {kr))F pq (k,t)dk 



(4.5) 



where Jo(kr) is same as (3.11) with n = 0. After substi- 
tuting (fOh in fl4.5Q and integrating we find 



/ r \ 2(1- D) 

S p (r,t)~(^) W) 2 ^ 



. 1 r Q t 

mthe ranges _ > _ > ^_ 

and 1 < £l Q t < Pe 1/3 . 



(4.6) 



V. STRUCTURE FUNCTIONS OF 
STATISTICALLY STATIONARY PASSIVE 
SCALAR 

To achieve a statistically stationary passive scalar field 
we may imagine that, as scalar patches take spiral shapes 
and decay, more scalar patches are introduced in to the 
flow as may well happen in an experimental setup in 
the laboratory. This procedure soon leads to a situa- 
tion where many scalar spirals coexist in the flow all in 
different stages of their evolution. Assuming the rate of 
regular injection of the scalar blobs to balance exactly the 
rate of scalar dissipation, we can expect to have a statis- 
tically stationary scalar field. In this case, the averaging 
over many spirals in different stages of their evolution 
(which is involved in the calculation of the generalised 
power spectra and structure functions) may be assumed, 
in the spirit of Lundgren JTcfl , to be equivalent to averag- 
ing over the life-time of a single spiral. Hence, to obtain 
the generalised power spectra of the statistically station- 
ary scalar we average the previous section's results over 
time in the range 1 < Qot < Pe 1 / 3 , which represents the 
life time of the spirals. The spiral structure lies in the 

range p < x n < L which implies flot < kL < yj 
This spiral range of wavenumbers together with the time 
range of the spiral gives the range of Slot over which we 
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can average for a given value of kL. This leads to a time 
averaged F pq (k,t) which takes the form 



F pq {k) ~ (kL)- 1 
where 1 < kL < Pe 1 / 3 , 

F pq (k) ~ (fcL)-(7-80) Pe 2(i-u) 
where Pe 1 / 3 < fc£ < Pe 1 / 2 . 



(5.1) 



(5.2) 



(5.1) and (5.2) are the result of averaging (4.4) over the 
time ranges 1 < Vl^t < kL and 1 < n^t < prjfyr respec- 
tively. These time ranges are determined by the respec- 
tive wavenumber ranges in (5.1) and (5.2). Finally (5.1) 
and (5.2) give the following structure functions 



S p (r) ~ constant + In ( 
where LPe" 1 / 3 < r < L, 



LPe- 1 / 3 



where 



6(1-13) 



LPe- 1/2 < r < LPe' 1 / 3 



(5.3) 



(5.4) 



I n (5.3|) and (5.4) S p (r) is a time average of S p (r,t) in 
(4.6). Note that C, p — with a logarithmic correction in 
the range LPe~ 1/3 < r < Land that ( p = 6(1-D) G]2,3[ 
in the dissipative range LPe- 1 / 2 < r < LPe- 1 / 3 . 

The structure functions S p (r) are not time-dependent 
and may be interpreted as characterising a scalar field 
in a statistically steady state achieved with an external 
large-scale source of scalar (scalar forcing). This scalar 
forcing may consist of regularly placing in the flow scalar 
blobs with large-scale smooth interfacial structure. In 
the spirit of Birkhoff 's Ergodic theorem |2^] we should 
expect Lundgren's time average assumption to be rele- 
vant for the calculation of structure functions when the 
scalar field is statistically stationary. 



VI. PHENOMENOLOGY 

In this section we extract the phenomenology under- 
lying the calculations and results of the previous sec- 
tions and show that the essential ingredient of this phe- 
nomenology are the locality of scalar inter-scale transfer 
(3.1) and the Lundgren time-averaging operation. In- 
deed, as we show in this section, (5.1) and (5.2) can be 
retrieved by a simple phenomenological argument based 
on these two ingredients. 

Let us return to the time-dependent wind-up of scalar 
spirals. As time proceeds, ie. as £lot — > Slo(* + 5t), then 
kL — ► kL + LSk because of the differential rotation(which 
amounts to local shear) in every steady vortex and the en- 
tire scalar spectrum is shifted towards higher wavenum- 
bers with time (see Figure 2). That is to say that the 



shearing advection to which the scalar patches are sub- 
jected in every steady vortex is such that the generalised 
power spectra obey 



F pq (kL + 5kL, fl t + n 5t)d(kL 
= F pq (kL,n t)d(kL) 



5kL) 



(6.1) 



3 



Q.„t+ Q.„8t 



kL 



FIG. 2. The time evolution of the scalar power spectrum 
F pq (kL,t) in a log-log plot. 

The amount of scalar variance in the wavenumber 
band d(kL) around wavenumber kL is simply transported 
to wavenumber band d(kL + 5kL) around wavenumber 
kL + SkL after an incremental time duration floSt (see 
Figure 2). As shown in p0|| , the distance I between con- 
secutive coils of the scalar spiral in every vortex at a 
distance r from the centre scales as 



I 



L 



1+1/ a 



(6.2) 



Letting time vary by a small amount St, the distance 
between two coils changes by 



81 
1 



St 

T 



the minus sign indicating that I is decreasing. Identify- 



ing k with ^f- for the purpose of equation (6.1) so that 
% = - y, it follows that (Oh becomes 



F„ 



kL I 1 + j 



dkL 



F pq {kL,n t)d{kL). 



The solution of this equation is 



F pq {kL,n a t) = L(kL)- x T pq 



ti t 
~kL 



(6.3) 



(6.4) 



where J- pq are arbitrary dimensionless functions. As in- 
dicated in figure 2 this form of the generalised spectra is 
valid in the limit Pe — > oo and in the wavenumber range 

fl t < kL < Jjfc and time ran S e 1 < ^t < Pe 1 / 3 . 
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Note that 1 < n t < Pe 1 / 3 < < Pe 1 / 2 . The 

inverse of fio* represents the decaying outer length-scale 



of the spiral range and the inverse of - 



represents the 



growing micro-scale of diffusive attrition. A wavenumbcr 
in the range 1 <C kL <C Pe 1 / 3 during the time-period 
1 <C flat <C Pe 1 / 3 does not have the time to be affected 
by diffusive attrition and only receives scalar variance 
activity from lower wavenumbers until Clot — kL. We 
therefore refer to 1 < fci < Pe 1 / 3 as the advective 
wavenumber range. The time averaged generalised power 
spectra in this range are given by 

F m {kL) = J d(n t)L(kL)-^ pq ^ 

(6.5) 

and T pq must be increasing functions of because the 
differential rotation's shearing process causes the power 
spectra to shift from small to l arge wavenumbers (see 
figure 2). Hence we retrieve (5.1), i.e. 



F pq {kL)~ (kL)- 1 

in the advective range 1 <C kL <C Pe 1 / 3 which is well 
defined in the limit Pe — > oo. 

In the advective-diffusive rang e Pe 1 / 3 < kL < Pe 1 / 2 
a wavenumber kL experiences the advection process from 

Clot = 1 until kL 



when molecular diffusion sets 

SZot 

in. Hence the time averaged generalised power spectra 
are given by 



Fpq(kL) = 



1 



Pe 



d(Sl Q t)L(kLY x T vq -4- 



D, t 



kL 



(6.6) 



in the advective diffusive range and using T pq 

\kL ) 



(see ([4.4[)) we retrieve (|5.2|), i.e 



F pq {kL) ~ Pe 2{1 - D) (kLY 



7+6D 



in the advective-diffusive range Pe 1 / 3 < kL < Pe 1 / 2 
which is well defined in the limit Pe — * oo. Note that the 
diffusive micro-length-scale LPe -1 / 2 is the Taylor mi- 
croscale of the scalar field (first introduced by Corrsin 
in 1951). Note also that 7 - 6D G]3,4[ and that the 
experimental results of || seem to show a steeper power- 
law wavenumber spectrum at wavenumbers larger than 
where the k~ l spectrum is observed. 



VII. CONCLUSIONS AND DISCUSSION 

In a two-dimensional isotropic and homogeneous col- 
lection of non-interacting compact and time-independent 



singular vortices with a large-wavenumber energy spec- 
trum E(k) ~ k~ a with 1 < a < 3, the structure functions 
of an advected and freely decaying scalar field have the 
following scaling behaviour in the limit where Pe — > oo 



, r vS(l-D) 

Sp(r,t) ~ [-n t) 



(7.1) 



~Pe~ ^ l ^ W~t and 1 < ^o^ < Pe 1/3 , and the 



fractal co-dimension P of the scalar interfaces is such 
that 1/2 < D < 2/3. 

By applying the Lundgren time-average assumption we 
obtain predictions for the structure functions of a statis- 
tically stationary scalar field in the same 2-D velocity 
field and the same limit Pe — > oo 



S p (r) ~ constant + In ( 



LPe- 1 / 3 



in the range LPe '' 3 C r « L and 
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(7.2) 



(7.3) 



in the range LPe -1 / 2 <C r <C LPe -1 / 3 . The logarith- 
mic term in (7.2) corresponds to k~ x generalised power 
spectra. It may be worth mentioning that the 2-D veloc- 
ity fields of Holzer and Siggia 0] where they observe a 
well-defined k~ x scalar power spectrum are replete with 
spiral scalar structures. 

Predictions of fc _1 scalar power spectra in the limit 
Pe — > oo have been made for scalar fields in smooth (ie. 
at least everywhere continuous and differentiable ) ho- 
mogeneous and isotropic random velocity field with arbi- 
trary dimensionality and time dependence by Chertkov 
et. al. p5fl w ho generalised and unified the results of 
Batchelor |26}] and Kraichnan (27). Experimental inves- 
tigations of the high Peclet number fc _1 scalar power 
spectrum have been inconclusive in 3-D turbulent flows 
even though Prasad and Sreenivasan have claimed such a 
spectrum in a 3-D wake [E8|. However fc _1 scalar power 
spectra have been observed at high Peclet numbers in 
numerical simulations of scalar fields in 2-D and 3-D 
chaotic flows p9|-|3l| and in 2-D velocity fields obeying 
the stochastically forced Euler equation restricted to a 
narrow band of small (integral scale) wavenumbers [B4| . 
More recently, k~ x scalar power spectra have been ob- 
served at Pe = 10 7 in 2-D or quasi 2-D statistically 
stationary turbulent flows in the same range where the 
velocity field's energy spectrum is k~ z by Jullicn ct al. 
M who have also observed logarithmic scalings in that 
range for all order structure functions (similarly to ([7.2]), 
but without the ability to establish the LPe -1 / 3 scaling 
factor and range). The theory of Chertkov et al. pa ] 
does not apply to this experiment because homogeneous 
and isotropic random velocity fields which are everywhere 
continuous and differentiable have energy spectra steeper 
than k~ A (see Appendix). The present paper's theory, 
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however, applies when the energy spectrum scales as k~ a 
with 1 < a < 3 but is limited to time-independent veloc- 
ity fields. Nevertheless the phenomenology developed in 
section VI also holds for time-dependent velocity fields 



and we now apply and generalise it to spatially smooth 
chaotic flows (and also, by the way, to frozen straining 
velocity field structures). 

The starting point of our phenomenology is the local- 
ity of transfer (^J). Pedrizzetti and Vassilicos ]2lJ have 
shown that inter-scale transfer in 2-D compact vortices 
is indeed local at a given scale when velocity gradients 
do not vary much in physical space over that scale. This 
is the case in the 2-D axisymmetric vortices considered 
in this paper but also in spatially smooth velocity fields. 
In a spatially smooth chaotic flow the distance / between 
successive folds of the scalar interface decays exponen- 
tially as determined by the largest positive Lyapunov 
exponent A, ie. l{t) ~ e~ At which implies y ~ — X5t. 
Applying the locality of transfer property we get 



F pq (k(l + X8t),t + St)dk(l + XSt) 



F pq (k, t)dk 



the solution of which is 



Fpq(k : t*) — k *Fpq 



(7.4) 



(7.5) 



This form of the generalised spectra is valid for Pe — > oo 
and as long as 1 < e xt < k, so that applying Lundgren's 
time-average operation from t = to t = A -1 lnfc gives 



F pq {k) = k- 1 



A 



Ink j 
[Hnfc]- 1 



T 

J j 



l»l 



dt 



(7.6) 



because F vq is an increasing function of £j— . Power spec- 
tra of scalar fields in spatially smooth chaotic flows are 
believed to scale as fc _1 in the limit Pe — ► oo |29-3lj but 
our theory predicts (fclnfc) -1 . This is a small correction 
to the spectrum but an exponentially large correction to 
the scalar variance in the limit Pe — > oo. 



VIII. ACKNOWLEDGMENTS 

MAIK and JCV acknowledge support from EPSRC 
grant GR/K50320 and from EC TMR Research network 
on intermittency in turbulent systems. MAIK also wishes 
to thank the Cambridge Commonwealth Trust, Wolfson 
College, Cambridge and DAMTP for financial support 
while this work was being completed. JCV acknowledges 
support from the Royal Society. 



IX. APPENDIX 

For a statistically homogeneous velocity field with ve- 
locity components Uj (x) we can define a correlation func- 
tion 



Rii(r) = («i(x)«j(x + r)) 
and its Fourier transform 

<M k ) = J dvR l3 {v)e-^. 

Incompressibility (V • u = 0) and statistical isotropy of a 
two-component or 2-D velocity field Uj(x), i — 1,2, imply 

kikj\ E(k) 



$y(k) = Sij - 



k 2 J irk 



where the average kinetic energy per unit mass of the ve- 
locity field is E = J °° dkE{k),k = |k|; E(k) is the energy 
spectrum of the velocity field. 

One dimensional energy spectra are defined as follows 

which, because of isotropy, arc completely characterised 
by a single component, say 4>ix{ki). 

From the above and a few standard manipulations one 
can get 



dx 



Vx 2 - 1 Ejxki) 
irx 2 0ii (fci) 



1 



and in a range of wavenumbers where both E(k) and 
0n (fci) are monotonically decreasing functions of k and 
k\ respectively, it follows that E(k) ~ k~ p if and only if 
0n(fci) - k^ p . 

The pivotal assumption in ]2aj is that the velocity field 
Mi(x) is Taylor-expandable up to at least first derivative 
terms everywhere in physical space. In particular, this 
means that iti(xi,0) is differentiable with respect to 
x\ everywhere on the x\ axis. If the first derivative of 
u\(x\, 0) with respect to x\ is also continuous everywhere 
along the x\ axis then the Fourier transform ui(fci) of 
u\ (x\ , 0) must decay faster than k^ 2 |^|| . If, however, the 
first derivative of ui(xi,0) is not everywhere continuous, 
then it is discontinuous either on a set of well-separated 
points or on a more pathological set of points which 
accumulate (and are therefore not well-separated) in a 
fractal-like or in a spiral- like manner (2^]3^] . In the case 
where discontinuities in the derivative field of iti(xi,0) 
are well-separated, u\(k\) decays as fcf 2 because the 
Fourier transform of well-separated discontinuities be- 
tween which the field is continuous decays as k^ 1 and 
the Fourier transform of the derivative of iii(xi,0) is 
equal to ikiUi(ki). In the other case where discontinu- 
ities are not well-separated and accumulate, the decay 
of u\(ki) can be anywhere between /cf 1 and fc-f 2 p4B3|, 
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but in this case U\(xx,0) cannot be considered to be 
differentiable at those points where discontinuities of its 
derivative accumulate. 

In conclusion, the differentiability of the velocity field 
everywhere in physical space implies that ui(fci) must 
decay at least as fcf 2 and therefore ^n(fci) ~ |ui(fci)| 2 ~ 
C(/cf 4 ) which in turn implies E{k) ~ 0(k~ 4 ). 

The condition E{k) ~ (D(k~ 3 ) stated in the conclu- 
sion of guarantees that the strain field is large-scale 
but not that the velocity field is differentiable. The spec- 
tral condition required to use the pivotal assumption of 
differentiability in || should in fact be E(k) ~ C(fc~ 4 ). 
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